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ABSTRACT 

Non-linear evolution of density fluctuations yields secondary anisotropies in the cosmic microwave 
background (CMB), which are correlated with the same density fluctuations that can be measured 
by weak lensing (WL) surveys. We study the CMB-WL correlation using analytical models as well 
as Y-body simulations. We show that an analytical model based upon the time derivative of non- 
linear matter power spectrum agrees with simulations. All-sky cosmic-variance limited CMB and 
WL surveys allow us to measure the correlation from non-linearity with high significance (50er) for 
'max = 10 4 , whereas the forthcoming missions such as Planck and LSST are expected to yield only 
marginal detections. The CMB-WL correlation is sensitive to the time derivative of structure growth. 
We study how this property may be used to constrain the nature of dark energy. While the expected 
constraints are not very strong, they may provide a cross check of results from other observations. 
Subject headings: cosmology: theory — large scale structure of universe — cosmic microwave back- 
ground 



1. INTRODUCTION 

Secondary temperature anisotropies of the cosmic mi- 
crowave background (CMB) provide invaluable informa- 
tion on the structure formation in the universe. The 
sources of anisotropies include the Sunyaev-Zel'dovich 
(SZ) effects by galaxy clusters as well as by reionization, 
the Rees-Sciama (RS) effect, and CMB lensing. 

The forthcoming Planck satellite and ground-based 
observ ations such as Atacama Cosmology Telescope 
CACT: lKosowskyl l2003ft and South-Pole Telescope (SPT; 
Ruh l et alJ 12004 are designed to measure temperature 
fluctuations at arc-minute scales, and thus are expected 
to detect some of the secondary anisotropies. 

The Rees-Sciama effect is, in principle, a unique probe 
of the ti me-variation of gravitational potential, as it is 
given by ijSachs fc Wolfdfl967t iRees fc Sciamal[l968h 



AT(n) 
T 



dr$'(hr, r), 



(1) 



where $ is the Newtonian potential, n is a unit direction 
vector, r is the conformal lookback time, and r* is r out to 
the photon decoupling epoch. The prime denotes d/dr, 
which is equal to —d/drj, where n is the conformal time. 

In the linear regime, <&' vanishes in a matter-dominated 
universe; however, when either curvature or dark energy 
dominates, $ decays in time and thus |$|' > 0. In the 
non-linear regime, on the other hand, <!> grows in time, 
|<f>|' < 0. Therefore, one expects $$' > on large scales 
where density fluctuations are still linear, and $$' < 
on small scales where fluctuations are non-linear. 

The auto-correlation of the RS effec t, both lin- 
ear and non-linear, has been studied (|Seliakl 119961 : 
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iTuluie et al.lH996l : lCooravll2002br ). The cross-correlation 
between the linear RS effect and large-scale struc- 
ture traced by galaxies has also been studied (e.g.. 



Crittenden fc Turoklll996l : lPeiris fc Sperge]||2000UCooravl 



2002at lAfshor di 2004), and detected for va rious galaxy 



surveys le.g.. iBoughn fc Crittenden! 12004 iNolta et all 
[200l lAfehordi et al.ll2004~ 

Much less attention has been given to the correlation 
of the non-linear RS effect and large-scale structure. As 
this effect changes the sign at the linear-to-non-linear 
transition scale, the signature of non-linearity is distinct. 

The RS effect is not the only thing that is correlated 
with the large-scale structure. The SZ effects, as well 
as radio and infrared point sources, also trace the large- 
scale structure. Fortunately, multi-frequency data enable 
us to separate these contributions that have unique and 
specific frequency dependence. 

In this paper, we calculate the cross-correlation of the 
RS effect and large-scale structure traced by weak lensing 
(WL) surveys. In particular, we simulate this directly 
using a Y-body code for the first time. 

We adopt the standard ACDM model with 
(q„n, n KQ , a s , h n,) = (0.26 0.74 ,0 .76, 0.7, 1), which is 
consistent with lSpergel et all (|2007| ). 



2. ANALYTICAL MODEL 

The amplitude of distortion in galaxy images due to 
WL is given by the so-called convergence, k, which is 
proportional to the projected density field along the line 
of sight. As the Newtonian potential, $, is negatively 
proportional to density field via the Poisson equation, 
the sign of CMB-WL correlation is given by — $$'; thus, 
negative on large scales and positive on small scales. 

The cross-correlation of the RS effect and n, C\ , is 



Q 



-IV 



drP^i (l/r; r) 



dz s p{z s )— r^, (2) 

r s r° 



where P$$/ (k; r) is the power spectrum of <&<!>' at r, z s 
is the redshift of source galaxies, and r s = r{z s ). Here, 
p(z)dz is the probability of finding galaxies between z 
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and z + dz. We use p(z) — Az 2 exp[— (z/zof], where A 
is a normalization factor determined by J Q p(z)dz = 1 
(|Efstathiou et al.lll991h . We consider two survey designs: 
(Model 1) Deep Survey, (j3,Zo) = (0.7,0.5), whose p(z) 
peaks at z ~ 2.2 with a broad distribution, and (Model 2) 
Shallow Survey, ((3, z ) — (2, 0.9), which peaks at z ~ 0.9 
with a narrower distribution. 
We calculate P$$/ from 

P **' (fc ' r) = 9 4a2°fc?° 4 [ P55 ' ( fc ' r ) " ^ p ss(k, r)] , (3) 

where = —dhxa/dr, and P55 and P55/ are the 
power spectrum of density fluctuations, S, and the cross- 
correlation of 8 and <5', respectively. 

How do we calculate P551I The continuity equation, 
5' = q where q = V ■ + 5)] is the momentum diver- 
gence, gives the exact relation, P$$> — Pg q . Here, Pg q is 
the density-momentum divergence cross spectrum. 

It might be tempting to use Pgg> = P$$/2, as the follow- 
ing might seem obvious: (Sk(r)5' p (r)) = \§^(5k(r)5 p (r)). 
However, this relation is exact only in the linear regime, 
for which the ensemble average (taken over initial pertur- 
bations) and d/dr commute. As they do not generally 
commute in the non-linear regime, we check whether this 
ansatz is a good approximation in the non-linear regime. 

We test this ansatz by comparing Pg q and Pgg/2. We 
use the 3rd-order perturbation theory (PT), as it allows 
us to study this delicate issue analytically. 

We expand Sk(r) in a series up to the 3rd order 
in the initial linear perturbation, <5i(fc), as 8k(r) — 
J2 n=1 D n (r)S n (k). We also expand the velocity di- 
vergence, = V • v, in Fourier space up to the 3rd 
order in the initial perturbation, 0i(k), as 0k(r) — 
^2 n -iD'(r)D n ~ 1 (r)6 n (k). Here, D is the linear growth 
factor, and S n and 9 n are of order 5™ and 6?™, respectively. 

With this expansion, we obtain 

P Sq (k,r) = D(r)D'(r)PU(k) 

+ D 3 (r)D'(r) [P™ (k) + P$ (k) + P^(k)] 

+B 1 g 21 g(k,p) + B 21 g 1 g(k,p)], (4) 

where Pgg(k) is the initial linear spectrum, Pgg(k) is the 

cross-correlation of 6i and 9j, and B% s (k,p) = (Si(k — 
p)9j(p)8i(k)) , whose explicit forms c an be calcula t ed in a 
straightforward manner by following lScoccimarrol (|2004T ). 

We t hen compare Ps a with t he 3rd order expression of 
P' ss /2 (|Bernardeau et alj|2002f ) 

\p l ss (k,r)=D(r)D'{r)PU(k) 

+ 2D*(r)D'(r)[P 22 (k) + 2Pti(k)]. (5) 

While equation ([4]) and (O are obviously identical to the 
first-order, it is not immediately clear whether Pg q = 
Pgg/2 holds in the non-linear regime. We find Pg q ~ 
Pgg/2 to the accuracy better than 10% for k < lh /Mpc, 
and 5% for k < 0.4/i/Mpc, at r = 0. 

In the right panel of FigureQ] we compare C; calculated 
from Ps q (dashed) and Pgg/2 (solid). We find that, in the 
non-linear regime (I > 1000) where the correlation turns 



positive, Ci from Pgg/2 is smaller by ~ 10%. This result 
suggests that we can use P551 = P' && /2 to compute P$$' 
from equation ([3]) reasonably accurately. 

However, our argument is still limited to PT, which 
breaks down for 5 > 1. In the same panel of Fig- 
ure [H the dotted lines show Ci from Pg S , where p^ is 
the fu lly evolved non-linear power spectrum (|Smith et al.l 
2003] ). The 3rd-order PT overestimates non-linearity i n 
Pss (jBernardeau et al.l 120021 : Ueong k Komatsul l2006t l. 
and thus PT predicts that the sign of Ci changes at 
smaller /. This result shows limitations of PT. 

3. ANALYTICAL MODEL VS iV-BODY SIMULATION 

To test P$$i « Pgg/2 in the non- linear regime, we com- 
pare the analytical model with iV-body simulations. 

We employ 512 3 dark matter particles in a volume of 
Lbox = 250ft.~ 1 Mpc and 40ft. _1 Mpc on a side. T he simu- 
lation s are performed by the GADGET-2 code ([Springe! 
2005). The initial conditions are generated by a stan- 
dard method using the Zel'dovich approximation. We 
dump outputs from z = 0.01 to 10 uniformly sampled in 
log(l + z). For each output redshift, we dump adjacent 
two outputs in order to calculate 5' and <&'. 

In the top left panel of Figured! we show the analytical 
model of Ci with Smith et al.'s Pgg (solid and long-dashed 
show Model 1 and 2, respectively) and the JV-body re- 
sults (open and filled symbols show negative and positive 
values, respectively). The agreement is good: the model 
describes the amplitude, shape, and cross-over at I ~ 800 
of Ci that are measured in the simulation. We therefore 
conclude that the ansatz, Pgg> w Pgg/2, is indeed accu- 
rate, up to I = 5000 where we can trust resolution of our 
iV-body simulation. 

How well are RS and WL correlated? In the bottom 
left panel of Figure [lj we show the 2-d correlation co- 
efficient, Ri — Ci/yJ CfCp s . We have used the same 
non- linear Pgg to calculate the power spectrum of conver- 
gence, Ci , while we have used the halo model approach 
(jCoorav k Shethl 120021 ) to calculate that of RS, Cf- S . 
The RS and k are strongly anti-correlated, Ri ~ — 1, 
in the linear regime. (It is not exactly —1 because 
of the projection effect.) The correlation weakens as 
non-linearity becomes important: Ri ~ 0.05 — 0.1 at 
10 3 < I < 10 4 for Deep (Model 1), and R t ~ 0.01 - 0.03 
for Shallow (Model 2). 

The weak correlation of C; is due to the fact that 
$ and $' are not correlated very well in the non- 
linear regime: the 3-d correlation coefficient, P3d(&) = 
P$$' (k) /-\/P# (k)P<s>' (k), reaches the maximal value, 
P 3D (0.5/Mpc) ~ 0.2, at z = 1, where the RS effect be- 
comes largest. The 2-d correlation, Ri, is even weaker 
than P3d(&) because of the projection effect and a mis- 
match between the redshift at which WL becomes largest 
(z ~ 0.5) and that for RS (z ~ 1). 

The weak correlation makes it challenging to measure 
the CMB-WL correlation from non-linearity. 

4. DETECTING NON-LINEAR REES-SCIAMA EFFECT 

Can we detect non-linearity, C/ > 0? The signal-to- 
noise ratio (S/N) of C\ is given by 

® <~ = /sky S 1 + (dr^/c 2 ) ' (6) 
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Fig. 1. — CMB-WL cross-correlation spectra, C;, for two WL survey designs: Deep (Model 1) and Shallow (Model 2). {Top left) 
The symbols show C; from the ./V-body simulation. The open and filled symbols show C; < and C; > 0, respectively. The solid and 
long-dashed lines show the fully non-linear model for Model 1 and 2, respectively, while the dotted lines show the linear theory predictions. 
Note that the linear theory predicts C; < at all I. {Bottom left) CMB-WL cross-correlation coefficients, Ri. {Top right) 3rd-order PT 
predictions. The solid and dashed lines show two ways of calculating C;, based upon Eq. J5]l and (|4jl, respectively. The fully non-linear 
calculations are also shown as the dotted lines (which are exactly the same as the solid and long-dashed lines on the top-left panel). {Bottom 
right) Difference between two ways of calculating C;. The fractional differences are at most ~ 10%. 



where / s k y is the fraction of sky observed, and Cf MB 

and Cf are the total (signal plus noise) power spectra of 
CMB (including primary, CMB lensing and RS) and WL, 
respectively. The noise spe ctra of CMB and WL survey s 
are given, respectively, by (|Knoxlll995t ISchneiderll2005l ) 

^yCMB = a 2.jf whm exp [i(i + l)0f whm / 8 In 2], (7) 

ATf = ^/n gall (8) 

where a plx is the temperature noise per pixel, #f w hm is the 
FWHM of a Gaussian beam, n ga i is the number density 
of galaxies observed in a WL survey, and <r 7 is the noise 
in shear measurements, including intrinsic ellipticities. 

We forecast S/N for the forthcoming surveys with the 
following parameters: (/ s ky, c 7 , n g /arcmin 2 )= (0.024, 
0.3, 20) for Canada- France-Hawaii Telescope Legacy Sur- 
vey (CFHTLS), (0.8, 0.1, 100) for Large Synoptic Survey 
Telescope (LSST), (/ sky , a pix , F whm) = (0.024, 4.4 M K, 
1'.7) for ACT, and (0.8, 1.7/2/4/3/xK, 10'.7/8'/5'.5) for 
Planck with 3 channels. 

In Figure^ we show the predicted S/N as a function 
of the maximum multipole, / max , for cosmic- variance lim- 
ited CMB and WL (Deep and Shallow) surveys, as well 
as for the forthcoming surveys: Planck correlated with 
LSST, and ACT correlated with CFHTLS. 

We find that S/N is totally dominated by the linear 
contribution (dotted lines) at I < 3000, and then becomes 
dominated by the non-linear contribution (dashed lines) 
at higher I, All-sky CMB and WL surveys can yield 
S/N ~ 50 (10) for Deep (Shallow) WL Survey, whereas 
1000 deg 2 surveys can only yield S/N ~ 7 (1). 

Once noise of the forthcoming surveys is included, 
however, S/N from the non-linear contribution be- 
comes small compared to the linear contribution. For 
Planck+LSST we find S/N ~ 1.5 for the non-linear, and 
6 for the linear. For ACT+CFHTLS we find S/N ~ 0.1 




lmax 

Fig. 2. — Signal-to-noise ratio {S/N) as a function of the max- 
imum multipole, l m ax- The dotted and dashed lines show the 
linear and non-linear contributions, respectively, while the solid 
lines show the total, for cosmic-variance limited CMB and WL sur- 
veys with all-sky (upper curves) and 1000-deg 2 sky coverage (lower 
curves). {Top) Deep WL Survey. The dashed-dotted lines show 
S/N expected from Planck's CMB data correlated with LSST's 
WL data. "PLANCK nl" shows the non-linear contribution only. 
{Bottom) Shallow WL Survey. Same as the top panel, but for 
ACT's CMB data correlated with CFHTLS's WL data. 

for the non-linear, and 0.7 for the linear. Therefore, we 
conclude that these forthcoming surveys are not expected 
to yield significant detection of non-linear RS effect. 

5. SENSITIVITY TO DARK ENERGY PARAMETERS 

As the linear RS effect vanishes during the 
matte r era, it is sensitiv e to d ark energy (DE) 
(e.g.. iBouehn fc Critteridenl l2004t iNolta et al.1 120041 : 
Afsh ordi et all 120041) . The non-linear effect measures 
the structure growth, which is also sensitive to DE 
(| Verde fc Sperge]||2002f) . The cross-over at / - 800, at 
which the linear and non-linear contributions cancel, is 
particularly a unique probe of DE. 
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Having evaluated the derivatives, we use the Fisher 
matrix analysis to calculate the expected constraints 
on DE parameters, p a — (£1\q, wq, wi), from the 
CMB-WL correlation. The Fisher matrix is given by 

F a p = X)i"2 X ( 5C, i/^a) Cov_1 ( 5C 'i/^/3) wn ere Cov is 
the covariance matrix (Eq. [5]) and Z max = 10 . 

The bottom panel of Figure [3] shows the expected joint 
constraints on wq and w\ from C; of all-sky cosmic- 
variance limited CMB and WL surveys. When Qao is 
not marginalized but fixed at the fiducial value, we find 
w = -1 ± 1.1 (0.4) and Wx = ± 4.1 (1.1) for Shal- 
low (Deep) WL Survey. The constraints are not very 
strong, unfortunately; however, the direction of degener- 
acy line, w\ ~ 2(wq + 1), is about twice as steep as that 
of a joint analysis of the baryon oscillations and Type la 
supernovae observations (|Seo fc Eise nstcin 20031). 



Fig. 3. — Sensitivity of the CMB-WL correlation to dark energy 
(DE). (Top) Derivatives of Ci with respect to the DE parameters. 
The solid, short-dashed, and long-dashed curves show dCi/dQ^Q, 
dCi/dwQ, and dCi/8w\, respectively. The thin dotted lines show 
those from linear theory. (Bottom) Fisher matrix forecasts for all- 
sky cosmic-variance limited CMB and WL surveys. The light and 
dark areas show the expected 1-<t constraints from Deep and Shal- 
low WL Surveys, respectively, with Qao marginalized. The solid 
and dashed lines show the same things with Oao fixed. 



The top panels of Figure [3] show sensitivity of Ci to 
DE parameters: dCi/dfl\o, dCi/dwo, and dC\jdw\. We 
use a simple form of DE equation of state, w(a) = 
Ph(a) I Pa(cl) , given by w(a) = wq + wx(l — a). The fidu- 
cial values are wo — — 1 and w\ = 0. 

We find dCi/dwo < and dCi/dwi < at all /. By 
increasing wq or wi, one makes w(a) less negative which, 
in turn, makes DE more important at earlier times. This 
does two things. On large scales where Ci < 0, it en- 
hances the linear RS effect, making Ci even more nega- 
tive. On small scales where Ci > 0, it reduces the growth 
of non-linear structure, making Ci less positive. In both 
cases we find negative derivatives. 

Dependence on Qao is more complex. While 
dCi/d ft ao < at most / can be explained by the same 
physics as above, we find dCij OHao > in the interme- 
diate I. This is due to the assumption of a flat universe: 
a larger f^Ao results in a smaller f2 m o, which alters the 
shape of Pss(k) by shifting it to larger scales. (A smaller 
!l m o delays matter-radiation equality.) This shift causes 
the intermediate / to behave differently. The behavior in 
the intermediate I has little to do with non-linearity, as 
we observe the same effect in the linear prediction. 



6. CONCLUSION 

We have studied the cross-correlation between CMB 
(the RS effect) and large-scale structure traced by WL. 
We have developed a simple analytical model based upon 
the time derivative of non-linear matter power spectrum, 
and tested its validity analytically with 3rd-order PT as 
well as numerically with A-body simulations. 

We have shown that all-sky cosmic-variance limited 
CMB and deep (shallow) WL surveys can yield a 50a 
(10c) detection of the non-linear CMB-WL correlation 
for Z max = 10 4 . The forthcoming surveys are not ex- 
pected to yield significant detections. We expect ~ 1.5a 
from Planck+LSST and O.lcr from ACT+CFHTLS. 

The change of the sign of Ci at the cross-over, I ~ 800, 
offers a unique probe of the nature of DE. Sensitivity 
of the CMB-WL correlation to DE turns out to be not 
very strong; however, as the direction of degeneracy on 
wo — wi is different from that of the baryon oscillations 
and Type la supernovae, it may provide an independent 
cross-check of the results from these observations. 
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